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ABSTRACT 



Aims. We use magnetic and non-magnetic 3D numerical simulations of solar granulation and G-band radiative diagnos- 
tics from the resulting models to analyse the generation of small-scale vortex motions in the solar photosphere. 
Methods. Radiative MHD simulations of magnetoconvection are used to produce photospheric models. Our starting 
point is a non-magnetic model of solar convection, where we introduce a uniform magnetic field and follow the evo- 
lution of the field in the simulated photosphere. We find two different types of photospheric vortices, and provide a 
link between the vorticity generation and the presence of the intergranular magnetic field. A detailed analysis of the 
vorticity equation, combined with the G-band radiative diagnostics, allows us to identify the sources and observational 
signatures of photospheric vorticity in the simulated photosphere. 

Results. Two different types of photospheric vorticity, magnetic and non-magnetic, are generated in the domain. Non- 
magnetic vortices are generated by the baroclinic motions of the plasma in the photosphere, while magnetic vortices 
are produced by the magnetic tension in the intergranular magnetic flux concentrations. The two types of vortices have 
different shapes. We find that the vorticity is generated more efficiently in the magnetised model. Simulated G-band 
images show a direct connection between magnetic vortices and rotary motions of photospheric bright points, and sug- 
gest that there may be a connection between the magnetic bright point rotation and small-scale swirl motions observed 
higher in the atmosphere. 
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1. Introduction 

One of the consequences of turbulent plasma movement in 
the solar photosphere is the appearance of horizontal vortex 
motions. These motions can be of significance for the gener- 
ation of magneto-hydrodynamic (MHD) waves whi ch prop- 
agate to the upper layer s of the solar atmosphere ([Parkeri 
119881 : iFedun et all 120091 : Uess et all 120091 ) . As a result of re- 
cent advancements in state-of-the-art instrumentation and 
\ observational techniques, it has now become possible to ob- 
serve small-scale vortices in the lower solar atmosphere. 

iBonet et al.l ([2008) showed the presence of vor- 
tex motions in photospheric G-band bright points 
with lifetimes comparable to those of the gran - 
ules. 



Wed emever-Bohm fc Rouppe van der Voortl (|2009[ ) ; 
■ IWedemever-Bohml (|2009l ) performed simultaneous G-band 
and Call 8542 A imaging to demonstrate the pres- 
ence of small- s cale swirl motions in the chromosphere. 
ICarlsson et al.l (|2010l ) also demonstrate a presence of the 
chromospheric swirls in their simulations, which include the 
chromospheric layer. 

A thorough investigatio n of the non-mag netic photo- 
spheric convection by IStein fc Nordlundl (| 19981 ) has shown 
the generation of vorticity by baroclinic fluid motions 
in the upper convection zone. Such motions, charac- 
terised by large and non-parallel gradients of density and 
pressure, occur near the edges of granules. Some stud- 
ies have also been performed on the vorticity generated 
by magnetic flux tubes rising to the solar surface from 
the deep sub-photosphere (e.g. lEmonet fc Moreno-Inse rtis, 
1998; E monet et alll200lh . The link between vortex gener- 



ation and magnetic field in th e upper photosphere was also 
noted by Vogl er et al.l (|2005f ). Generation of Alfvenic vor- 
tices by the interaction of compres sible plasma wi t h field - 
aligned o bstacles was stud i ed by iGruszecki et al.l (|2010l ). 
Recently, iKitiashvili et al.l (|2010f ) used numerical simula- 
tions to demonstrate the significance of vortex motions and 
vortex dragging for the creation of pore-like magnetised 
structures in the photosphere. 

Here we present a more detailed study of the photo- 
spheric vorticity. We focus on its origins and connection 
to the photospheric magnetic field and granulation dynam- 
ics. We show a direct correspondence between the photo- 
spheric vortices that correspond to G-band bright point 
motions, with the strong intergranular magnetic field. We 
also demonstrate a presence of vortex motions in the upper 
atmosphere, which may be connected to the chromospheric 
swirls. 

In Section 2 we provide a brief description of the numer- 
ical model, simulation setup and the results obtained. The 
output of the simulations is analysed in terms of the vortic- 
ity equation in Section 3. Section 4 describes the radiative 
diagnostics and observational consequences of vortex mo- 
tions in the photosphere, while in Section 5 we summarise 
our conclusions. 



2. Simulations 



We use the MuRAM code (| Voder et all I2005D to per- 
form the simulations, which has been s uccessfully used 
for a wide range of solar a pplica t ions ([S chiissler et al 



tor a wide range ot solar applications (bcnussler e t al. 
200$ IShelvag et all 12004 120071 : ICheung et alT l^OOS 
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Fig. 1. Geometry of the domain. Stream lines (red) are plotted over the three-dimensional structure of the vertical com- 
pone nt of the magnetic field (blue-green). The image is generated with the VAPOR 3D visualisation package ([Clvne et all 
120071) . 



Pietaril aGraham et all 120091: iRempel et all |2009; 
lYelles Chaouche et all 120091 : iDanilovic et all l2010h . The 
code solves large-eddy radiative three-dimensional MHD 
equations on a Cartesian grid, and employs a fourth-order 
central difference scheme to calculate spatial derivatives. A 
fourth-order Runge-Kutta scheme is used to advance the 
numerical solution in time. Hyper diffusivity sources are 
used to stabilise the solution against numerical instabilities 
and to account for physical processes which are not resolved 
by the numerical grid. Non-grey radiative energy transport 
is included in the code using short-characteristics and 
opacity binning techniques. A non-ideal equation-of-state, 
taking into account the first-stage ionisation of the 11 most 
abundant elements in the Sun, is also included. 

The numerical domain has a physical size of 12 x 12 Mm 2 
in the horizontal direction, 1.4 Mm in the vertical direc- 
tion (?/), and is resolved by 480 x 480 and 100 grid cells, 
respectively. The lower boundary is transparent, allowing 
the plasma to travel in and out of the domain. The upper 
boundary of the domain is closed, however, it allows the 
horizontal motions of the plasma and the magnetic field 
lines. The side boundaries of the domain are set to be pe- 
riodic. The domain is positioned in such a way that the 
visible solar surfac^B is located approximately 600 km be- 
low the upper boundary. 

Our starting point for the simulations is a well- 
developed non-magnetic (B = 0) snapshot of photospheric 
convection taken at t ~ 2000 s (about 8 convective turnover 



1 In this paper we refer to the visible solar surface as a hori- 
zontal geometrical layer which is physically close to the optical 
layer of radiation formation. We assume that the geometrical 
properties of the analysed features do not significantly change 
within this layer. 



timescales) from the initial plane-parallel model. A uni- 
form vertical magnetic field of B y = 200 G has been in- 
troduced at this stage, and a sequence of 147 snapshots 
recorded, containing physical parameters of the model, such 
as velocity and magnetic field vectors, temperature, density, 
pressure and internal energy. The sequence covers approx- 
imately 40 minutes of physical time, corresponding to ~5- 
10 granular lifetimes. During the simulation, the magnetic 
field is advected into the intergranular lanes by the con- 
vective plasma motions, and the maximum field strength 
rises from its initial value of 200 G to a few kilogauss 
in the intergranular lanes. The sequence we obtained is 
long enough for the relaxation of the model, since the time 
needed for the magnetic field redistribution is about .1-0.2 
hours ([Voder et all 120051 : I Voder fc SchiisslerL I2007D . 

A three-dimensional rendering of the velocity and mag- 
netic field is shown in Fig. [TJ Stream lines (red curves) 
reveal a large amount of vortex motions in the upper pho- 
tosphere. These vortices appear primarily in the intergranu- 
lar lanes and coincide with regions of strong magnetic fields 
and downflows. 

The modulus of the horizontal velocity components is 
shown in Fig. [2j where the top panels correspond to a level 
close to the visible solar surface level, and the bottom ones 
to a height in the domain close to the temperature min- 
imum. Images on the left correspond to the initial non- 
magnetic model, and those on the right to the fully devel- 
oped magnetic snapshot. It is evident from the figure that 
small-scale vortex structures have formed in the magne- 
tised model in the upper photosphere (bottom-right panel). 
These structures are not seen in the non-magnetic model 
nor at the visible solar surface. The contours in the bottom- 
right panel of Fig. [2j which bound the granular regions 
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Fig. 2. Horizontal cuts of the modulus of the horizontal components of velocity in the domain. The cuts are taken 
approximately at the visible solar surface (top panels) and in the upper photosphere (bottom panels). The left panels 
correspond to the initial non-magnetic snapshot, while the right panels correspond to the well-developed magnetic model. 
Vortex motions are clearly visible in the bottom-right panel. The contours on the bottom-right panel bound the regions 
where the vertical component of magnetic field B y < 30 G. 



where the vertical component of magnetic field B y < 30 G, 
clearly demonstrate that the vortices in the upper photo- 
sphere are co-spatial with the magnetic field concentrations 
in the intergranular network. 



3. The vorticity equation 

defined 



as the curl 
Similar to 



of 



velocity. 
Stein fc Nordlundl 



Vorticity is 

uj = V „ . 

(11998ft: iMoreno-Insertis fc Emonetl dl996); 

lEmonet fc Moreno-Insertisl ([1998ft : lEmonet et all (j2001[ j. 
we write the vorticity equation as the curl of the momentum 
equation of the MHD system: 



D uj 

p__ = ("-V)v. 



V- X Vp q - 
P 



+V x 



-J x B 



P 



(1) 



By substituting the current vector J with V x B and 
combining the magnetic pressure term with the gas pres- 
sure, we obtain the following equation: 



Dt p 



D uj 1 

P—+ — = (<*>' V)V- V- X V(p g +Pm) + 

+Vx -(B- V)B, 
P 



(2) 



where D/Dt represents the full (material) derivative, uj 
is the vorticity vector, p and p g are the plasma density and 
pressure, v is the velocity vector, B is the magnetic field 
vector, and p m = B 2 /2 is the magnetic pressure. Note that 
the magnetic field is normalised by a factor for conve- 
nience. Viscous dissipation of vorticity has been neglected 
in this equation. The material derivative in the left side of 
Eq. (pQ) is expressed in terms of Euler derivatives as: 

D uj duj 

^7 = ^ + w(v - v) + (v - v)w ' (3) 

Eq. (|2]) can be rewritten by separating the magnetic 
and non-magnetic terms and decomposing the last mag- 
netic term into two: 
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Fig. 3. Vertical component of vorticity uj v at the visible solar surface (upper panels) and in the upper photosphere (lower 
panels). Non-magnetic (left) and magnetic (right) snapshots are shown. 
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The right-hand side of Eq. (|4]) shows the different phys- 
ical mechanisms associated with the generation of vorticity. 
T\ is the vortex tilting term, T2 is responsible for the hydro- 
dynamic baroclinic vorticity generation, T3 represents the 
magnetic baroclinic vorticity, and T4 corresponds to vor- 
ticity generated by the magnetic tension. The quantity in 
square brackets within T3 represents the deviation of the 
magnetic field configuration from a potential field. 

Horizontal cuts of the vertical component of vorticity 
are shown in Fig. [3l The upper panels show the vorticity 
maps at a level close to the visible solar surface, while the 
lower panels correspond to the upper photosphere. Non- 
magnetic and magnetic snapshots are shown in the left and 
right panels, respectively. 

In the non-magnetic case, the vorticity is generated by 
the hydrodynamic baroclinic term and is mostly randomly 



directed without any internal structure of the vortices in 
the intergranular lanes. However, we notice a change in the 
structure of the vortices after we incorporated the magnetic 
field into the numerical box. Once the magnetic field has 
been advected into the intergranular network, both posi- 
tive and negative polarities of the vorticity coexist within a 
vortex, with both clockwise and counterclockwise motions 
appearing. An increase in the amplitude of the vertical com- 
ponent of vorticity by a factor of about 5 is also observed. 

In order to determine the origin of the vorticity, the 
terms in Eq. (|4|) need to be analysed separately for the 
whole simulation sequence. The vertical component of vor- 
ticity, which corresponds to the horizontal vortex motions, 
is of primary interest. Thus, a mean value of the modulus 
of the vertical vorticity component is a good measure of the 
amount of vertical vorticity generated in the model. 

Fig. [H shows the evolution of the mean of the modu- 
lus of the vertical vorticity component at different heights 
in the domain, together with the time dependences of the 
^-components of T\ through T4 terms of Eq. (|4|). The left 
panels correspond to the upper photosphere (500 km above 
the visible solar surface), the middle panels the approxi- 
mate level of the visible solar surface, and the right panels 
correspond to the convection zone (650 km below the visible 
solar surface). 
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Fig. 4. Dependence of < \u y \ > (top plots) and of different terms of the vorticity equation (bottom plots) on time at 
different heights in the domain. The first column corresponds to the upper photosphere (500 km above the approximate 
visible solar surface level), the second column is the visible surface level, and the third column corresponds to the 
convection zone (650 km below the visible surface level). 



An inspection of the left panels of Fig. H] reveals that 
the vorticity in the upper layers of the model is produced 
by the magnetic field, rising from its non-magnetic value 
of 0.0025 s _1 to 0.01 s _1 within the first three minutes of 
the simulation. During this phase, the magnetic field gets 
almost completely transported in to the intergranular lanes 
by convective motions of plasma ([Vogler et al.L l2005f ) . The 
amount of vorticity produced in the photosphere (middle- 
top plot) remains roughly the same, experiencing some ini- 
tial decrease, which may be connected to the suppression 
of plasma motions by the initially uniform magnetic field, 
before it increases again at t = 2 — 4 min. The behaviour 
of the vorticity at the bottom of the domain is opposite to 
what is observed at the top: in the first 5 minutes of the 
simulation, the amount of vorticity has decreased by almost 
a factor of 2. 

The evolution of vorticity shown in Fig. H] can be ex- 
plained by analysing the relative importance of the vertical 
components of the T\ through T4 terms in Eq. (|4]). In the 
upper photosphere, the term which corresponds to the vor- 
ticity generation by the magnetic tension (T4, solid black 
line) experiences a sharp rise, and after the initial phase 
of the simulation (4 min) takes the largest value among 
the other terms. The same behaviour is observed for the 
magnetic baroclinic term T3 (red dashed line), with an am- 
plitude which is a factor of 3 smaller than that of T4. The 
hydrodynamic baroclinic term T2 (blue dash-dotted line) is 
a factor of 2 smaller than T3. Thus, the baroclinic motions 
of the fluid (in the hydrodynamical sense) do not make a 
significant impact on the vortex generation in the upper 
photosphere. 

A different picture is observed at the photospheric level. 
Here, the amount of vorticity produced by the hydrody- 
namic baroclinic motions of the fluid is similar to the 



amount of vorticity, generated by the magnetic tension. In 
the initial stage of the simulation, T4 and T2 behave in an 
opposite way. This is caused by the processes of magnetic 
field redistribution. Initially, the magnetic field is uniform, 
and suppresses the plasma motions both in the granules 
and in the intergranular lanes, thus decreasing the vortic- 
ity generation by hydrodynamic baroclinic term T2, while 
T4 experiences a sharp increase due to the formation of the 
magnetic flux concentrations. After the initial stage, when 
the convection pushes the magnetic field out of the gran- 
ules, both the magnetic- type vortices and hydrodynamic 
vortices can be produced. 

In the convection zone, (right panels in Fig. H]), where 
the plasma /3 is high, the hydrodynamic baroclinic term T2 
is the primary source of vorticity. Although T2 decreases 
significantly from its non-magnetic value within the first 
few minutes of the simulation, it remains a factor of 2 larger 
than the magnetic tension term T4. 

We note that T4 plays a significant role for the vorticity 
generation both in the convection zone, where the plasma 
f3 is large, and in the upper photosphere, where plasma /3 
is small in the magnetic field concentrations. This term has 
no "hydrodynamic" equivalent and represents a physical 
mechanism for vorticity generation which is separate in na- 
ture from the conventional baroclinic vorticity generation 
processes. The magnetic baroclinic term T3 does not show 
any significant influence on the generation of vorticity. The 
significance of the hydrodynamic baroclinic vorticity gen- 
eration term increases with depth. 

3 mHz acoustic oscillations are also present in the do- 
main. Both the vertical component of vorticity uo y and T4 
show signs of these oscillations (see upper left plot of Fig. |4j) 
after the simulation has passed its initial stage. This find- 
ing, together with the connection of the photospheric vor- 
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tices to the magnetic field, confirms the idea that the oscil- 
lations leak through magnetic concentrations to the upper 
photosphere. An observational search for oscillatory signals 
in vorticity may be possible. 

4. Radiative diagnostics 

Photospheric bright points corre spond to regio ns of 

strong intergra n ular magneti c fields ([ Schussle r et all [2003; 
ICarlsson et all l2QQ4t IShelvag et all I2004D and may be 
subject to vortex motions. Recent high spatial resolution 
observations indica te that this may indeed be the case 
([Bonet et al.l ([20081 )). The analysis we present includes the 
radiative di agnostics in the G-band. Using the method s de- 
scribed by iFrutigerl (|200Q[ ). iBerdvudna et al.l ([2003D and 
IShelvag et al.l (|2004l ) , we have computed G-band images for 
all 147 sequential simulation snapshots. The direct effect of 
the magnetic field on the absorption line profiles has not 
been included in the calculations. The sequence of images 
allowed us to not only to find and track the vortex mo- 
tions of the photospheric G-band bright points, but also 
to study their appearance during the initial stages of the 
simulation while the magnetic field was being advected to 
the intergranular lanes. An inspection of the images made 
possible to identify vortex motions associated with mag- 
netic bright points (MBPs) in the simulated photosphere. 
An example is shown in Fig. [5] and confirms the connec- 
tion between the photospheric vortex motions, rotation of 
MBPs and vorticity in the upper layers of the simulated 
photosphere. Three consecutive simulated G-band images 
of an MBP (top panels), together with vorticity maps at 
a height corresponding to the visible solar surface (mid- 
dle) and in the upper photosphere (bottom) are shown 
in the figure. All images are taken at the same horizon- 
tal position. The evolutiona ry track of the M BP is very 
similar to the observations of lBonet et al.l ([2008). The bot- 
tom i mages clearly resemble the evolution of chromospheri c 
swirl ( We demever-Bohm fc Rouppe van der Voortl l2009f h 
despite the data being obtained somewhat lower than the 
level of Call core formation. This fact suggests that the 
chromospheric swirls may be connected to the vortices, gen- 
erated in the solar magnetic photosphere, however, a fur- 
ther investigation is needed to provide a rigorous proof for 
that. 



5. Conclusions 

In this paper, we have examined the processes leading to 
vorticity generation in the simulated magnetised photo- 
sphere. We have shown that large amount of vorticity in 
the photosphere is formed as a result of the photospheric 
plasma interaction with the magnetic field in the inter- 
granular lanes. The amount of vorticity generated due to 
the magnetic field is a factor of 4 larger than that gener- 
ated by baroclinic motions in convecting photospheric non- 
magnetic plasma. By the appropriate decomposition of the 
magneto-hydrodynamic vorticity equation we defined two 
vorticity equation terms, which are connected to the mag- 
netic field: the first resembles the baroclinic term in hydro- 
dynamics, while the second contains the magnetic tension 
and does not have a direct hydrodynamic equivalent. We 
have demonstrated that it is only the latter term which 
is mostly responsible for the generation of vorticity in the 



upper photosphere. Conversely, the importance of the baro- 
clinic hydrodynamic term increases with increasing geomet- 
rical depth. We also note the appearance of oscillations in 
the magnetic tension term in the upper photosphere. These 
may be signatures of the 3 mHz lower-photospheric oscilla- 
tions which leak to the upper layers through the intergran- 
ular magnetic field concentrations. 

Using radiative diagnostics with the G band, we con- 
firmed that MBPs are subject to rotary motions in the in- 
tergranular lanes and are magnetically connected to the 
vortices in the upper photosphere. 

The large volume occupied by the intergranular vortices 
suggests their significance for the energy balance of the solar 
atmosphere. Further analysis is needed on the connection 
of the photospheric vortices with the chromospheric swirls, 
with the upper chromosphere and corona regions, and on 
torsional wave excitation in these regions of the solar at- 
mosphere. 

Follow-up investigations will focus on the details of the 
physical mechanism, which leads to the creation of neg- 
ative and positive vorticity signs in a magnetic vortex. 
Extending the simulation box to include the chromosphere 
with a full non-LTE t reatment of the radiative diagnostics 
([Carlsson et all l2010f h will also be the subject of a future 
investigation. 
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Fig. 5. Simulated G-band images of an MBP (top panels), vorticity at the visible solar surface (middle panels), and the 
vorticity in the upper photosphere (bottom panels). 
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